o 



Charmonium from Classical Pure SU(3) 
Yang-Mills Configurations 

O. Oliveira* R. A. Coimbra''" 
O ■ Centre de Fisica Computacional, Departamento de Fisica, 

^ ' Universidade de Coimbra, 3004-516 Coimbra, Portugal 

February 2, 2008 

(N 

CN ' Abstract 

' A generalized Faddeev-Niemi ansatz for the gluon field is discussed. In 

, its simplest parametrization, the ansatz allows a solution of the classical 

• SU(3) Yang-Mills equations. From these solutions a confining potential 

' for heavy quarkonia is defined. The investigation of charmonium spectra 

. proves that the potential is able to reproduce the experimental spectra at 

' the level of 3%. Moreover, charmonium leptonic and strong decays are 

I investigated. The results are in fair agreement with experimental figures 

fH , and are in line with other quark model calculations. 



1 Introduction and motivation 



Quantum Chromodynamics (QCD) is the theory describing the interaction be- 
tween quarks and gluons. QCD being a non-abelian gauge theory, besides the 
^ , quark-gluon interaction includes gluon-gluon interactions. It is expected that 

^ ' these interactions are responsable for the usual baryons and mesons together 

with exotic states like, for example, glueballs and gluelumps. 

In QCD, the basic building blocks are the quark and gluon fields. However, 
it appears that nature has a preference for baryons and meson states, i.e. colour 
singlet states, rather than single isolated quarks or exotic states. Quarks seems 
to be permantly confined within hadrons. The theory provides hints about the 
quark and gluon confinement mechanism, but a fully explanation is still lacking. 

The simplest picture suggesting quark confinement is the singlet potential 
computed from Wilson loops [T1|3I21E]- This potential can be viewed as the 
energy between two static colour charges. For large quark separations, the 
potential grows linearly making it impossible to separate two quarks. In this 
sense, the singlet potential is a confinement mechanism applicable to heavy 
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quarks. For light sector, even the idea of describing the quark-quark interaction 
via a non-relativistic potential is questionable. 

Besides the Wilson singlet potential, two popular confinement mechanisms 
are the dual Meissner effect O El El El and the vortex condensation picture 
Pl IIUI im 1121 [T^ . In both pictures, confinement is explained invoking special 
types of gluon configurations, namely the condensation of magnetic monopoles 
in the dual picture and the percolation of colour vortices in the condensation 
scenario. Although, a considerable amount of work has been done in favour of 
both mechanisms, their precise confirmation is still missing. 

The above described scenarios do not exhaust possible explanations for quark 
confinement. Many other mechanisms can be found in the literature. Typically, 
these mechanisms look for certain types of gluon configuration which, in princi- 
ple, explain why we are not able to detect a free quark. Implicitly, they assume 
that the origin of confinement is within the pure gauge sector. The success of 
the quenched approximation in lattice QCD further supports such an idea. 

The interaction between heavy quarks, i.e. charm and bottom, can be stud- 
ied with non-relativistic quark models Hopefully, the interquark potential 
^1 El El El El should be derived from QCD. In practice, the potentials are 
either phenomenological or, at best, QCD-motivated. The idea of defining an 
interquark potential is appealing and, if one could solve the theory, the potential 
should be related to a given gluon configuration as the Coulomb potential in 
QED. However, so far, the known solutions of the classical field SU(3) gauge 
theory, i.e. instantons, do not suggest quark or gluon confinement. 

In this paper we construct solutions of the classical field SU(3) Yang-Mills 
theory, in Minkowsky spacetime, which lead to confining nonrelativistic poten- 
tials^. The classical solutions are obtained after the introduction of a gener- 
alized Cho-Faddeev-Niemi-Shabanov ansatz which, in Landau gauge, allows a 
parametrization of the gluon field in terms of two vector fields and a scalar field. 
For such a gluonic field, the classical field equations are essentially a system of 
coupled QED like equations for the two vector fields and a massless Klein- 
Gordon equation for the scalar field. The solutions of the coupled equations can 
be divided into QED like solutions, i.e. plane wave solutions, and exponential 
growing fields. The later ones will allow to identify a confining potential for 
heavy quark systems. Although the potential is given by a multipole expansion, 
in this work only the I = term will be considered. In order to compute its 
parameters, the potential is compared with the singlet potential computed from 
the Wilson loop. It comes that, in the range 0.2 - 1 fm, the new potential fol- 
lows closely the Wilson loop behaviour, with a maximum deviation of about 50 
MeV. The differences being both the small (0(l/r)) and large (©(exp (Ar)/r)) 
distance behaviours. Moreover, the potential is characterized by a single mass 
scale A = 228 MeV, computed together with the other parameters by fitting the 
singlet potential in the region 0.2 - 1 fm. Being close to the singlet potential, the 
charmonium is well described by the new potential. Moreover, assuming a pure 
vector interaction one can define a spin dependent potential which, together 

1 See 1201111] for related work. 
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with central potential, is able to reproduce experimental spectra with an error 
of less than 3%. The leptonic and sam hadronic decays of charmonium are also 
investigated. The model is is fair agreement with the experimental figures. 

The paper is organized as follows, in section 2, for completeness, the gen- 
eralized Cho-Faddeev-Niemi-Shabanov ansatz proposed previously is discussed, 
giving particular attention to its simplest parametrization. In section 3, a po- 
tential for heavy quarks is motivated and compared with the singlet potential 
computed from the Wilson loop. The potential parameters are computed fitting 
the singlet potential in the region 0.2 - 1 fm. In section 4, the charmonium spec- 
tra is calculated with the new potential, including the effect of spin-dependent 
forces. The recently discovered charmonium states and their possible quantum 
numbers assignment are discussed. Furthermore, discussion of the hadronic 
decays and leptonic decays of the 1 is performed. Finally, we draw the con- 
clusions and plans for future work in the last section. 



2 Classical Gluon Fields 

For SU(3) the lagrangian density^ reads 

^^_lj,a^^a^. (1) 

where 

f;, = d^Al - d^Al - gfabcAlAt , (2) 
and A° are the gluon fields. The classical equations of motion, 

d^F^" + ig [A^ , F^"-] = 0, (3) 

being a set of non-linear partial differential equations are quite difhcult to solve 
and to compute a solution of it is usual to introduce an ansatz. 

Following a procedure suggest in I23| , let us consider a real covariant constant 
field n". From its definition, it follows 

D^n'^ = d^^n'^ + ig{F-)abAln'' - 0, (4) 

where f are the generators of representation to which n belongs. If one wants 
to parametrize the gluon field in terms of n, the scalar field should have, at 
least, as many color components as A^. The adjoint representation fullfils such 
condition and it will be assumed that n belongs to this representation, then 

fbac , (5) 
d^^n- + gfabcn'Al = 0. (6) 

Multiplying the last equation by n'^ it follows 

a^Kn") =0. (7) 

^In this work we use the notation of )22| . 
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From now on it will be assumed that n°n° = 1. 
Defining the color projected field 

- n-Al (8) 

one writes the gluon field as 

= n'^C^ + X;, (9) 

where X is orthogonal to n in the sense 

n^X^ = 0. (10) 

In order to establish the gauge transformation properties of the various fields, 
we consider infinitesimal gauge transformations, 

U{x) ^ e-'"(^) = l-iu){x) (11) 
5AI = U^uj'^ + fabcOj'Al (12) 

Sn"- = fabc^^'n'. (13) 

For Cp, the transformation rule looks like a "projected abelian" gauge transfor- 
mation 

5C^ = in'^a^c^". (14) 

For X°j^, the transformation property is non-linear in the field n. It follows from 
®, (ini and ^ 

<5X^ = i {S""" - n'^r?) df^Lu" + fabcCj'X'^ (15) 

and 

S{n^X^^) = 0. (16) 

In order to build an usefuU ansatz, one can go back to equation (O, replace 
the gluon field by the decomposition lO, multiply by faden'^ and try to solve 
for the X field. After some algebra we get 

2 

fabcn'df.n" - gff^^" ^ g {dabeddce - ddbedace) Xl^n'^n'^ = (17) 



after using the SU(3) relation 

2 
3 



where 

^abc 



iTr(A'^{A\A^}) (19) 
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and A" are the SU(3) Gell-Mann matrices. Equation p7|l allows us to write 

K = ^fabcn%n^ + y; (20) 

and 

n'^X^ = n'^Y/^ = 0. (21) 

The first term of is a generalized Faddeev-Niemi ansatz The simplest 
non-trivial n field with norm one can be parametrized as 

n" = S^^sine + S^^cosO. (22) 

Then, with the above definitions equation H17|l implies the following relations 
between the Y fields 

Y^ = Y^ cote, (23) 

- ^^d,9, (24) 

= Y^ = Y^ - = 0, (25) 

which together with the color space orthogonality condition (|21|l gives the gluon 
field 

= n'^C^ - ^-^d^e + S-^'B^- (26) 

in the last equation we defined = Y^. The corresponding gluon field tensor 
components are 

= ^"C^- + S^^B^,,, (27) 

where 

C^_iu = d^Cn — 9i/C^ , Bfj_i, = df^Bi, — d^B^ . (28) 

For the above gluonic configuration, the classical action is a functional of the 
vector fields and _B^, 

C = - i (C2 + 6^) , (29) 
and the classical equations of motion (PJ become QED-like equations 

a^c^, = 0, = 0. (30) 

Moreover, the hamiltonian density and the spin tensor are given by the sum of 
the contributions of two abelian-like theories associated with and B^ fields. 
Note that the pure gauge theory at the classical level, is independent of 9. How- 
ever, the inclusion of fcrmions implies a coupling 0-fermions. The coupling with 
9 is associated with the first three Gell-Mann matrices and, for configurations 
with = B^ = requires only ; no coupling with the third color component 
and the coupling to the first two color components have opposite signs. 
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The computation of classical solutions of the equations for the gluon config- 
uration considered requires gauge fixing. For the Landau gauge^, the full set of 
equations is 

dt'Cf, = 0, d^Bf, = 0, (31) 
d^eCf, ^ 0, d^d^e = 0. (32) 

These together with the classical field equations have plane wave solutions char- 
acterized by a four momenta k such that — 0, the fields C^t and i?^; are 
polarized perpendicularly to fc, 

C^{k,X) = ec(fc,A)e~''=^ k^ec{k,X) ^ 0, (33) 

B^{k,X) = eB(A:,A)e-^'=^ k''eB{k,X} = 0, (34) 

0{k) = Ooe-''"'' , (35) 

and the third color component being polarized along k. Our aim is to identify 
possible solutions of the classical field equations (|31(l . (|32|l . H3()|l which can sug- 
gest confining solutions. We will not proceed with the discussion of the solutions 
of the above set of coupled equations but, instead, discuss a particular class of 
solutions. 

For the pure gauge theory, the condition for a finite action/energy solution 
does not constraint 9. Let us consider the simplest non-trivial configuration one 
could think of, namely = = 0. This particular configuration has zero 
energy and the field equations are reduced to a massless Klein-Gordon equation 
for 9 which can be solved by the usual separation of variables. Indeed, writing 
9{t,r) = T{t)V{r), it comes 

_ (36) 
T{t) ~ V{r} - ^ ' ^"^^^ 

where A is mass scale which is invariant under rescaling of the gluon field. The 
solutions with < are the usual free field solutions. For A = 0, the gluon 
field is linear in time and the spatial part is the usual solution of the Laplace 
equation. If A^ > 0, 

Tit) = Are"^* + Btc-^' (37) 

and writing 

V{r) = J2 Vi{r)Yi„,{n) (38) 

where Yi^iQ) are the spherical harmonics, the functions Vi{r) can be writen 
in terms of the modified spheri 
Ki+i/2iz) kind, where z ~ A?-, 



in terms of the modified spherical Bessel functions of the 1^^ h+\l2{z) and 



Vl (r) = Aim 1= \- Birn 1= . (39) 



■^For the Coulomb gauge, the first two equations are shghtly more comphcated but in the 
third and fourth equation one should replace 9 by V, together with other obvious modifica- 
tions. 
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Figure 1: Minima of as function of A. 



These solutions have the same scale for the temporal and spatial parts and, 
asymptotically, are real exponential functions. For example, the lowest multi- 
pole solution is 

, , , sinh(Ar) e"^'' , , 

Voir) = A ^ + B 40 

r r 

and the associated gluon field is 

Al = A(e^*-6Te-^*) Vo{r) , (41) 
P = _(eA* + 6Te-^*) Vyo(r). (42) 



3 A Non-relativistic Potential for Heavy Quarko- 
nium 

The coupling of the classical field configuration discussed in the previous section 
to the fermion fields requires only the Gell-Mann matrix, i.e. 

T'^ya^^V'; (43) 

ip should be understood as a quark Dirac spinor. Assuming that it makes 
sense to describe the heavy fermion interaction via a potential derived from 
the above gluon configuration, from the structure of one naively expects a 
bound state, let us say the first color component, an unbound state, the second 
color component and a free particle solution, the third color component. Of 
course, the usual quark model picture for hadrons is recovered only after gauge 
transforming all fields. 

The classical solution discussed in the last part of the previous section has 
the same spatial and time scale A = 228MeV, i.e. its spatial scale is L = 0.865 
fm, while its temporal scale T = 2.886 x 10~^^ s. 



7 



For charmoninum, typical time scales ranges from r = 1/91 KcV ^ = 7.2 x 
10-21 s for J/ip{lS) up to T = 8.4 x lO'^^ s for ^(4160). 

For bottomonium the time scale ranges from r ~ 1.5 x lO^^^ s for T(15') up 
to r = 6.0 X 10-24 s for T(10860). 

Typical time scales for charmonium and bottomonium are much larger than 
the time scale of the classical configuration. If quarks propagate in the back- 
ground of the above classical solution, the effects of propagation should be small, 
i.e. the quark interactions are almost contact interaction in time. Therefore, in 
first approximation, the interaction can be parametrized as 



6{t' - t)V{x' - x)ij{t', x')i!{t, x) . 



(44) 



On the other end, typical wave lengths of charmonium and bottomonium 



mesons are 1/3000 MeV" 



0.07 fm and - 1/10000 MeV" 



0.02 fm. 



respectively, and they are much smaller than L ^ 0.9 fm. From the point of view 
of the classical configuration, the mesons are point like structures, i.e. they can 
be described by a point like wave function ^(r, t). Furthermore, if one assumes 
that charmonium and bottomonium are non-relativistic quark systems, then it 
seems reasonable to describe the mesons via a gauged Schrodinger equation with 
a non-relativistic potential given by H38|l . In first approximation one can forget 
the contribution from the vector potential of the classical solution. 

The non-relativistic potential H38() is not spherical symmetric. However, if 
only the lowest multipole multipole H40() is considered, one is back to the central 
potential picture. From now on, we will assume that the confining potential for 
heavy quarkonium is given by the lowest multipole of the spatial part of and 
that the Schrodinger equation for heavy mesons is, in first approximation, 



(45) 



where (/) is a heavy meson field and m the reduced mass of the quarkonium 
system. 

The potential ipO]! is coulombic for small quark distantes 



r 



and confining for large distances 



A 



Voir) ^ ^e-; 



(46) 



(47) 



since the potential grows exponentially with distance, it can be viewed as a soft 
bag model. For large inter-quark distances, the wave function goes to zero as 




Ar 

T 



for r ^ 1 , 



(48) 



with the spatial extension of the quark field becoming smaller for heavier quarks. 
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Figure 2: = Voir) - Vs,ngiet{r). 



The first step towards a proper definition of Vq as a potential for heavy quark 
systems is the computation of its various parameters, namely A, B and A. This 
can be done by minimizing the square difference between the new potential and 
the singlet potential computed from Wilson loops. For the singlet potential we 
used the results from |3J at /? = 6.4. The 



=0.2/m 



1 2 



Voir) - Vsingletir) 



(49) 



as function of A can be seen in figure ^ The curve shows an absolute minimum 
and we take the potential parameters as their values at this point 

A = 11.2542671, B = - 0.70113530875, A = 228.026 MeV (50) 

when r is measured in MeV^^. The difference between Vo and the lattice po- 
tential is reported in figure |21 

It is curious that the value of the only energy scale in the solution is close to 
standard values for Aqcd- We have no interpretation for this result. Probably, 
it is connected to the way the parameters are defined, i.e. to the fit of Vo to the 
lattice singlet potential. 



4 Charmonium - the spectra 

The potential Vb(r) is a first approximation to the full interaction. From the var- 
ious possible corrections, it will be assumed that the spin-dependent part gives 
the dominant contribution. Assuming a pure quark-quark vector interaction, 
the spin potential (see, for example, or ^3]) is given by 

Vspvn_^ Vso+Vten + Vss, (51) 

^Assuming that the short distance behaviour of the inter-quark potential is dominated 
by the one-gluon exchange contribution, then B = — 4as/3. The strong coupling constant 
associated to the fitted value for B is = 0.526. 
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where the spin-orbit Vgo, tensor Vten, and spin-spin Vss potentials are 

Vso = :j^-^L-S, (52) 

2mg r dr 



Vss = :i VVos- -s,, (54) 

where L is the orbital angular momentum, sq the spin of the heavy quark, toq 
its mass and 



Sfpn — 4 



3 (sQ • f) {sq -f) - Sq 



(55) 



is the tensor operator. With the above definitions, the potential verifies the 
Gromes constraint that arises from the boost invariance of the QCD. 

In the following, the mesons states are classified according to the usual spec- 
troscopic notation n^^~^^Lj, where S is the total spin and J the total angular 
momenta. The Schrodinger equation is solved for Vo(r) with the wave function 
^nim('r) = Rniif) Yijn{0, (/)) , whcrc n IS the principal quantum number, I and m 
are the orbital angular momentum and its projection, Rni(r) is the radial wave 
function and Yim{6,(j)) are the spherical harmonic functions. The total wave 
function is then built by coupling "^nim{r) with the spin wave function. The 
contribution of the spin-dependent potentials will be included perturbatively. 

We further define 

i.e. the radial wave function and its derivatives at the origin. The various 
R^jiO) are required to compute production rates and decays rates. 



4.1 Quark mass 

The charm quark mass, rric, was tunned to reproduce the Xc2i^P) — J/'4'{^S) 
mass difference. Using the mesons mass values from the particle data book j2()| . 
the mass diference is 

M (xc2(l-P) ) - M (j/1'(lS') ) = 459.34 ±0.11 MeV (57) 

and rric = 1870 MeV. We have studied other possible definitions for rric but 
although rric changes, the meson spectra is essentially the same discussed here. 



4.2 The spectra 

The charmonium radial wave function and its derivatives, for various states, are 
reported and compared with other quark model calculations in tabled 

For S wave mesons and for Vq, the wave function at the origin is essen- 
tially constant, except for the state IS". For the other potentials and for the 
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QCD(BT) 


Power law 


Logarithmic 


Cornell 


IS 


1.739 


0.810 


0.999 


0.815 


1.454 


2S 


1.261 


0.529 


0.559 


0.418 


0.927 


3S 


1.283 


0.455 


0.410 


0.286 


0.791 


2P 


0.248 


0.075 


0.125 


0.078 


0.131 


3P 


0.394 


0.102 


0.131 


0.076 


0.186 


3D 


0.222 


0.015 


0.026 


0.012 


0.031 



Table 1: Charmonium radial wave function and its derivatives at the origin. 
All quantities are in GeV^'+"^. The first column reports the value computed 
using Vo{r) and all others are taken from The references for the various 

potentials are: QCD(BT) 16 , Power law [HI, Logarithmic jH], Cornell 



states reported in tabled ^(0) decreases when the principal quantum number 
increases. This is a major difference between Vq and the other potentials. For 
Vq this means that the decays widths which are sensible to the wave function 
at the origin, such as the leptonic decays, should scale approximately with the 
inverse of meson mass squared. For the other potentials this scaling behaviour 
is not verified. Looking at i?(0), the potential which is closer to Vq is the Cornell 
potential. 

For P mesons, the first derivative of the wave function at the origin in- 
creases with the principal quantum number for the QCD inspired potentials Vq , 
QCD(BT) and Cornell. The largest increase reported in table^happens for Vq. 
For the other two potentials, |i?*^^-'(0)p is essentially constant. 

For the only D meson reported in table ^ Vq shows a second derivative 
at the origin an order of magnitude larger than all other potentials. Thus, Vq 
enhances, for example, the leptonic decay of the D meson states when compared 
with all other models. 

The full theoretical charmonium spectra is reported in table EJ together 
with experimental spectra, including the recently discovered charmonium states 
X{3872) 1^1^, Z(3930) (xc2(2P) in the latest version of the Particle Data 
Book 1^) ISO], X(3940) ini and r(4260) 32 . The experimental values are 
from the Particle Data Book with the exception of the new states. The 
meson masses are defined according to 

M(meson) = 2mc + Enr + ^cc , (58) 

where Enh is the eigenvalue of the Schrodinger equation and Acc = —3223 
MeV is a shift introduce to reproduce the J/^E" experimental mass value. This 
mass shift, Acc, can be viewed as a non-perturbative contribution which is not 
accessible within the model. 
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On overall, one can find good agreement between the experimental numbers 
and the model predictions. The deviations reported in table |21 are below 3 % 
level. 

The quantum numbers of various charmonium states have never been mea- 
sured experimentally. In the following, we comment on our results and possible 
quantum number assignments. 

4.2.1 V(3836) 

The experimental information on ■0(3836) is scarce and this state no longer ap- 
pears in the last particle data book j^. Only two experiences |i34j have 
seen signs for this meson. The measure of the quantum numbers of ^'(3836) was 
discussed in |34|. Based on the comparisation with theoretical mass predictions, 
the authors suggest that it should be a state. Following the same reason- 
ing, from the values of the mesons mass and considering the already excluded 
quantum numbers, this state can be identified with a l^D^. The mass difference 
between the prediction and the experimental figures is 5 MeV (0.1%). 

4.2.2 X(3872) 

This particle was discovered by Belle collaboration |^ in the decay B~ 
X{3872)K-, with X J/V'TT+TT", and confirmed by CDF DO ^ and 
BABAR collaborations. Experimentally, only the charge parity C = -I- is 
well established with the data favouring |4l)[ 1411 the assignment J^^ = 1++. 
A J^^ — 2++ is not ruled out but seems to be unlikely. In our calculation, l^D-z, 
2^Pi, 1^F2 mesons have masses close to the mass of X{3872) and with C = +. 
The assignement with 2^Pi state shows a deviation from the experimental mass 
of 67 MeV (1.7%). 

4.2.3 The 3940 MeV mass states 

Belle collaboration reported on the three possible new states with masses around 
3940 MeV: A:(3940) 31 , y(3940) ^31 and Z(3930)/xc2(2P) EOI with measured 
masses of 3943 ± 6 ± 6 MeV, 3940 ±11 MeV and 3931 ± 4 MeV, respectively. 
The ^(3940) was observed in the process e+e~ J/tpX, y(3940) was seen 
as a resonance in the decay B KttttttJ /ip and Z(3930) was observed in the 
reaction 77 D^D , D~ . For the Z particle, the helicity distribution of the 
final states particles favors a J = 2 assignment. In the update version of Particle 
Data Book "5^^ it appears as a Xc2(2P) state with mass value 3929 ± 5 ± 2. 

Looking at the mass values computed with the model under discussion, the 
natural candidates are 2P and 1^F2 levels. For Z, assuming that it has J = 2 
and C = -|-, it should be a 1^F2 charmonium state. The theoretical computed 
mass shows a deviation from the experimental number of 3 MeV (0.0%). 

In what concerns the quantum numbers, for the remaining two states there is 
no experimental information. Looking at the mass values, the only candidates 
for X(3940) and r(3940) are 2^Fo and 2^Pi states. Since the identification 
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of X with a Xci(2P) state has problems [5^], one assumes that it is a 2^ Pi 
state. The theoretical prediction is 111 MeV (deviation of 2.8%) lower than 
the experimental value. For y(3940) the mass suggests a Xco(2P) state. Such 
state has positive parity state, in agreement with the prediction of |S7]. Note 
that according to these authors, X(3940) should be one of the following states 
Xci(2P), hc{2P) or r]c{3S). Their analysis favors the 3^5*0 state. For Y the 
theoretical prediction is 16 MeV (deviation of 0.4 %) above the value quoted in 
the particle data book. 

4.2.4 y(4260) 

This broad resonance was observed in e^e~ 7/5_R7r"''7r~ J/t/j with a mass 
of 4.26 GeV .32 , meaning that it has the quantum numbers of the photon 
J^^ = 1 . From table [51 the closest mass state is the 3'^ Si, which has a mass 
96 MeV (deviation of 2.3 %) lower than the measured value. 

4.2.5 V(4040), V(4160), ^(4415) 

The states -0(4040), V'(4160), -0(4415) were observed in e+e^ — > hadrons 45 . In 
what concerns their quantum numbers, there is no experimental information but 
it is usual to assume that they are vector particles. In what concerns theoretical 
predicitions, there aren't 1~~ states around the mass value of ^"(4040) and 
\['(4415). With quantum numbers 1^^ we can assign the ^1/(4160) with 2^Di. 

From table 121 and looking only at the mass values, there a number of states 
which can be associated with ^^(4040) and ^'(4415). For the first state, 2^P2 and 
l^i^3 are good candidates and for ^'(4415) the candidates are 3'^ Pi and 3^Pi. 
Note that, given the difference in mass to the lowest 1 meson, 942 MeV for 
\P(4040) and 1324 MeV for ^(4415), a possible interpretation for both states 
being that they are gluonic excitations rather than quark excitation states. In 
this sense, *(4040) should be a bound state of J/* and the J^^ = 0++ /o(980), 
while ^-(4415) a bound state of J/* and the J^^ = 0++ /o(1324). Given the 
lack of experimental information, hopefully, only after studying their production 
and decays will be possible to identify the associated quantum numbers. 

4.3 Comparisation With Other Quark Models and With 
Lattice QCD 

In order to compare our spectra with other charmonium calculations, we repro- 
duce table II from '56" including our figures in tableland updating the particle 
data book figures. In table (31 BGS are the results of a simple non-relativistic 
quark model GI is an updated version of the "relativised" Godfrey-Isgur 
model EFG is the relativistic quark model of Ebert, Faustov and Galkin 
|54| . Cornell are the results of the original Cornell model including effects to 
the coupling to open charm virtual states ^Hj- CP-PACS [3Si and Chen 
are two lattice calculations in the quenched approximation. For comparisation, 
in the first column we report the Particle Data Book values I^HI- Note that in 
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table, the particle identification assumes the "Standard" notation, i.e. 3^Si is 
identified with ^-(4040), A^Si with ^-(4415) and Z(3930) with Xc2(2P). From 
the table it is clear that Vq as well as the lattice simulations underestimate the 
hiperfine splittings — ^ 5*0. For the lattice this is believed to be due to the 
quenched approximation - see, for example, [44) . 

Following 56 , in table^we report the average model errors using the values 
reported in table |3| Note that, the particle identification in table 13 does not 
favors our calculation. Nevertheless, the model predictions are in line with other 
quark models and provides an overall agreement with the experimental spectra 
similar to lattice calculations. Of course, an analysis of the meson decays is 
required to confirm or not the model as a valid model for charmonium. 



5 Charmonium Decays 

In order to further test the heavy quark potential under discussion, we have 
computed leptonic and hadronic decay widths. 

5.1 Leptonic Decays 

For the leptonic decays we follow the van Royen-Weisskopf approach and 
assume that QCD corrections can be included perturbatively. Then, the com- 
putation of the decay widths requires only the knowledge of the wave function 
at the origin and its derivatives .47^ . 

For n^Si states, the leptonic width Fg+g- is given by 

^fl-i'Si) - ^\Rnsm' (59) 

where Mqq is the mass of the "^Si state, e<j is the electric quark charge (in units 
of |e|), a the fine structure constant and Rns{0) the radial wave function of the 
bound state at the origin. 

For n^Di states, the leptonic width is |60[ I61| 

e'.-m) - ^l«'^'(0)P (60) 

(2) 

where Rlj^^ stands for the second derivative of the radial wave function at the 
origin. 

In table |S1 we report our predictions for the leptonic widths. In order to 
avoid the problem of the estimation of QCD corrections, the theoretical widths 
were computed relative to the J/"^ width. The figures reported in table |S1 use 
the particle data book value for 

r(j/q,^g+e-) = 5.55 ± 0.14 ± 0.02 keV. (61) 
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In table our predictions for n^Si widths are compared with other models. 
All values where computed relatively to the J/^' width and using th experimen- 
tal value for the meson masses as reported in the Particle Data Book [201 • All 
models give reasonable estimates for the ^'(25') leptonic width but for the re- 
maining states the leptonic width are clearly overestimated, with Vq providing 
always the largest figures. For Vb, 5* (4040) does not fit into the spectra and 
therefore the number reported should not compare well with the experimental 
figure. For \1/(4160), at least in some models, it is expected a large mixing with 
an S-state which could explain the differences between the theoretical prediction 
and the experimental measure. 



5.1.1 Mixing Effects in Charmonium 

For J^^ = mesons, the theoretical predictions of the leptonic widths com- 
puted using the heavy quark potential Vq are reported in tableEl For '0(25') the 
model prediction is larger than the experimental number by a factor of 15%. 
The model overestimates the widths of -0(3770) and 0(4160) by a large factor. 
The deviation can be explained, at least partially, in terms of mixing between 
states. In the quark model under discussion, the tensor interaction is respons- 
able for the mixing of the S-wave and D-wave central potential states. Indeed, 
■0(25) a 2^Si state, 0(3770) a l^Di state and r(4260) identified with a 3^51 
state, \['(4160) identified with a have the same spin and the same total 

spin. Moreover, the theoretical mass predictions for such states differ by 29 
MeV and 9 MeV, respectively; see table El 

The perturbative calculation of the mixing induced by the spin-dependent 
potential improves only very sligthly the agreement between theory and ex- 
periment. Therefore, the differences are not due to the mixing between the 
single channel states considered here but mainly to the coupling to open charm 
channels. In this article we will not perform a coupled channel analysis. In- 
stead, a phenomenologically estimate of the mixing for these J^*^ — 1 states, 
parametrizing the mixing via a single parameter |51| 

1^(3770)) = cos^ll^Di) + sin6'|235'i) , (62) 
1^(25')) = -sin6'|l3i:>i) + cose* |23S'i) , (63) 



The leptonic widths are given by 



2„2 



(^(^^™)) = M^[^(3770)] 



5 



sin6li?2s(0) + — cos6ii?f^(0) 
2v2m^ 



2 



and 

re+e-(^(25)) 



4a2e2 



M2[V'(25)] 



cos0i?2s(O) - smOR^^liO) 



(64) 



(65) 



where a is the fine-structure constant and Cc the charm electric charge. Similar 
expressions hold for y(4260) and ^'(4160). If the QCD corrections are identical 
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for the two states it cames 



M2[7^(3770)] r,+e- (V'(3770)) 



1.123 sin 61 + 0.089 cos 



1.123 cos6' - 0.089sin^ 
0.092 ±0.012 



Fitting these ratios to the experimental values yields the solutions 9 
1.0°, —21.4° ± 1.0° and the leptonic decay widths become 



(66) 
(67) 

12.4° ± 





no mixing 


with mixing 


Exp. 


^(23) 
V'(3770) 


2.84 
0.02 


2.62 
0.23 


2.48 ±0.06 

219+° °2^ 



with all widths in keV. Clearly, the mixing improves the theoretical estimations. 

Assuming that y(4260) is a S^Si state which mix with ^-(4160) a 2^Di 
state, or reversing the assignments, and adjusting the mixing angle to repro- 
duce the experimental width Fg+e- (^'(4160)) — 0.83 ± 0.07 keV one obtains 
Fg+e- (y(4260)) — 1.42 keV in both cases. This calculation is resumed in the 
following table. 



state 


assignment 




Te+e- 








no mixing 


with mixing 


Exp. 


3'^ Si 


*(4160) 


2.28 


0.83 


0.83 ±0.7 


2^Di 


r(4260) 


0.04 


1.42 




3^51 


r(4260) 


2.16 


1.42 




2^Di 


^'(4160) 


0.05 


0.83 


0.83 ±0.7 



5.2 Hadronic widths 

In what concern hadronic decays, for S'-wave mesons one can consider the fol- 
lowing processes^ 

n^So — > •yj, gg. (68) 
The theoretical estimates of widths, including QCD corrections, are given by 
(see, for example 03) 

F( 5o ^ 77) = 3e,a (^1 - j (69) 

Spor n^Si mesons, one can computed theoretical estimates for partial widths for decays 
into 777, ggg and ggy. However, typically, ^Si states couple to a large number of hadronic 
decay channels, which makes the comparisation of the theoretical estimates with experimental 
figures quite difficult. Moreover, the QCD corrections are large and probably the factorisation 
used in the theoretical calculations does not hold. Indeed, typically, the theoretical figures are 
quite far away from the experimental values. For these reasons, we do not discuss hadronic 
widths for '^Si mesons. 
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with the replacemente nric — M {meson) /2. 

For u^Sq states, again considering only ratios to the ground state meson, 
and using the following experimental values for T{ric 77) — 6.7jig'g keV and 
r('7c gg) = Ttot = 25.5 ± 3.4 MeV we obtain the widths reported in 
table [7| For j]c{2S) the theoretical estimated total width agrees well with the 
measured value. 

For P wave mesons the decays we consider the following decays 

n^P2 — > 77, gg, 

n^Pi — > qqg , 
n^Pa — > 77, 55 

n^Pi — > 77, ggg, ggi- 

(71) 

The expressions for the decay widths relevant to P states, without QCD cor- 
rections, are 

n'Po -> gg) = 6a^ \^^E^ (72) 
rfPo-77) = 27e^a2^^5iiM (73) 
r(3Fi ^ qqg) = ' ^ ^ ^' ln(mQ < r >) (74) 



5 mj 



n'P2 --gg) - ^ (75) 



r( P2 77) = ^e,a —-4 (76) 

r(^^'i -> .9.97) = f e^^r(iPi^.g55) (78) 
5 ^ as 

where (r) — drru'^i{r). The expressions for nP states are from except 
H78() which is from |^. The predicted widths are reported in table |S| In the 
calculation we have as = 0.29 ± 0.03 (see [521), neglect QCD corrections and 
used M{state)/2 rather than ruq in the formulas, except in the logarithm, where 
we used AI {state) < r >. Note that these formulas are rough estimates of the 
partial widths for the corresponding annihilations processes. 

The rates for annihilations of the '^Dj states into three gluons (via color- 
singlet operators) are given by jBSJ |^ 

2 

lOa^ ^d(O) 

r{^Dj ^ ggg) = —^Cj — ln(4mQ < r >) (79) 
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where Cj = f, 1, 4, for J = 1, 2, 3. 

The two-gluon annihilation rate of the state is given by |52| : 



:3 R';,iO) 



2 



2al 



,2 a' 



n'D2 ^ gg) 



3 




(80) 



The values in table |H1 are similar to the ones calculated by except for the 
state l^Di where they used M{l^Di) ~ 3872 MeV, instead of the experimental 
value. 

6 Results and Conclusions 

In this paper it is proposed a generalized Faddeev-Niemi ansatz for the gluon 
field - see equations (0 , l|20|l and (|21|l . The ansatz writes the gluon field in terms 
of two vector fields and a real scalar field in the adjoint representation of SU(3). 
In it simplest parameterization, A^^ is given by (|26|l and the associated classical 
equations of motion become abelian-like equations. Moreover, the hamiltonian 
and spin tensor are the sum of two abclian-like contributions from the two vector 
fields. In Landau gauge, the solutions of the classical equations of motion can be 
writen as plane waves, with a time-like four momenta, and with the having 
longitudinal and transverse polarizations. A particular class of solutions with 
zero classical energy is built. These configurations are asymptoticaly growing 
exponential functions in space and in time and are caracterized by a unique 
mass scale. From these solutions, a potential for heavy quark phenomenology is 
suggested and investigated. In the range of 0.2 - 1 fm, the new potential follows 
essentially the singlet potential but differ from it for small and large inter-quark 
separations. The minimization of the square difference between the potentials 
gives A = 228.0 MeV. 

The potential, including perturbatively the spin dependent contributions, 
is then used to investigate charmonium. In what concerns the mass spectra, 
the description is precise up to 3% level. Typically, the difference between the 
theoretical estimation and measured mass is well below de 100 MeV. However, of 
the known charmonium states, two mesons do not fit in the model prediction, 
namely ^'(4040) and ^'(4415). Possible interpretations of these particles as 
charmonium states is suggested but not discussed in detail. A full study of the 
possibilities rised here for these particles is under way. Despite of the success 
of the spectra, the hiperfine splitting cames two low in the model. This is 
seen also in lattice simulations and, in our case, probably cames from fitting 
the new potential to the lattice single potential. We are currently trying to fit 
the potential parameters to charmonium spectra. This program is very time 
consuming but we hope to be able to report the results in a reasonable time 
scale. 

A first study of the mesons decays was performed. In what concerns S'-wave 
leptonic decays, the new potential describes well the widths of the lower states. 
For excited states, the identification is not clear and the comparisation with 
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experimental numbers requires further investigations. The leptonic decays show 
that, for certain states, mixing is crucial. We have computed perturbatively the 
mixing induced by the spin-dependent potential. However, the mixing induced 
by the spin-dependent potential only changes slightly our results. On one hand 
this shows that our calculation is consistent. On the other hand shows that a 
coupled channel analysis is required to reproduce the observed behaviour. We 
are currently starting such a program. 

In what concerns the hadronic decays, the model predictions are in reason- 
able agreement with total widths estimates for ^So and P-wave mesons. 

The results for the charmonium spectrum and decays suggest that the po- 
tential and, hopefully, the classical configurations discussed here can teach us 
something about heavy quarkonium. At least, they seem to provide a starting 
point for a more elaborated calculation, namely a coupled channel calculation, 
including possible exotic contributions. We are currently involved in exploring 
such possibilities and extend the use of the potential to other heavy quarkonium. 
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Table 2: Charmonium spectrum up to 5000 McV. The table includes only L < 4 
states. The quantum states with * where used to set the model parameters (rric, 



State 


PDG 


Vo 


BGS 


GI 


EFG 


Cornell 


CP-PACS 


Chen 


J/*(1^5i) 


3096.916 ±0.011 


3097 


3090 


3098 


3096 


3095 


3085 ± 1 


3084 ± 4 


2980.4 ± 1.2 


3075 


2982 


2975 


2979 


3095 


3013 ± 1 


3014 ±4 


1''(2'\5'i) 


3686.093 ± 0.034 


3659 


3672 


3676 


3686 


3684 


3777 ± 40 


3780 ± 43 




3G38 ± 4 


3()33 


3()3() 


3()23 


3588 


3084 


3739 ± 46 


3707 ± 20 




4039 ± 1 


4164 


4072 


4100 


4088 


4225 


- 


- 




4135 


4043 


4064 


3991 


4110 


— 


— 


*(4«S'i) 


4421 ± 4 


4669 


4406 


4450 


— 


4625 


— 


— 




4638 


4384 


4425 




4460 




— 


X2(1^P2) 


3556.20 ±0.09 


3556 


3556 


3550 


3556 


3523 


3503 ± 24 


3488 ± 11 




3510.66 ±0.07 


3462 


3505 


3510 


3510 


3517 


3472 ± 9 


3462 ± 15 


Xo(l'Po) 

ftc(l'Pl) 


3414.76 ±0.35 


3372 


3424 


3445 


3424 


3522 


3408 


3413 ± 10 


3525.93 ±0.27 


3478 


3516 


3517 


3526 


3519 


3474 ± 40 


3474 ± 20 


X2(2^P2) 


3929 ± 5 ± 2 


4048 


3972 


3979 


3972 


— 


4030 ± 180 


— 


XI (2^ A) 




3938 


3925 


3953 


3929 


— 


4067 ± 105 


4010 ±70 


Xo(23Po) 




3832 


3852 


3916 


3854 


— 


4008 ± 122 


4080 ± 75 


/ic(2iPi) 




3959 


3934 


3956 


3945 


_ 


4053 ± 95 


3886 ± 92 


X2(3^P2) 




4546 


4317 


4337 


— 


— 


— 




Xi(3^Pi) 




4422 


4271 


4317 


— 


— 


— 




Xo(33Po) 
/ic(3iPi) 




4300 


4202 


4292 


— 


— 


— 






4446 


4279 


4318 




- 


— 








3841 


3806 


3849 


3815 


3810 


— 


3822 ± 25 


*2(1^£'2) 




3755 


3800 


3838 


3811 


3810 


— 


3704 ± 33 


r/c2(liD2) 


3771.1 ±2.4 


3688 


3785 


3819 


3798 


3755 


— 


— 




3754 


3799 


3837 


3811 


3810 


— 


3763 ± 22 


*3(2^i?3) 

*2(23i?2) 




4327 


4167 


4217 


— 


4190 


— 


— 




4230 


4158 


4208 


— 


4190 


— 


— 


*i(23Di) 


4153 ±3 


4155 


4142 


4194 


_ 


4230 


_ 


_ 




4230 


4158 


4208 


_ 


4190 


_ 


_ 


X4(1^P4) 




4114 


4021 


4095 


- 


- 


- 


- 


X3(1'P3) 




4012 


4029 


4097 








4222 ± 140 


X2(1'P2) 




3932 


4029 


4092 










/ic3(l'P^) 




4007 


4026 


4094 








4224 ± 74 


X4(2^P4) 




4603 


4348 


4425 










X3(23P3) 
X2(23P2) 
/lc3(2lP3) 




4492 


4352 


4426 












4406 


4351 


4422 












4487 


4350 


4424 










*5(1'G5) 




4382 


4214 


4312 










«'4(1''G4) 




4260 


4228 


4320 










«'3(1''G3) 




4161 


4237 


4323 














4253 


4225 


4317 











Table 3: Quark models and lattice charmonium spectra. 
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BGS 


GI 


EFG 


Cornell 


CP-PACS 


Chen 


1.7% 
0.6% 


0.3% 
0.1% 


0.6% 
0.2% 


0.4% 
0.2% 


1.8% 
0.7% 


1.5% 
0.6% 


1.4% 
0.5% 



Table 4: Average Model Errors. The first line is the average error. The second 
line is the square root of the sum of the relative deviations squared, divided by 
the number of states. 





4?(o) 


2 

^theo 


assign. 


r 

^ exp 


2^51 
3^51 

43^1 


1.261 
1.283 

1.350 


2.84 
2.26 
2.28 
2.16 
1.90 


*(25) 

*(4160) 
*(4260) 


2.48 ±0.06 
0.83 ±0.07 


3^Di 


0.031 
0.101 
0.222 


0.02 
0.05 
0.08 


^'(3770) 
*(4160) 


2iq+°°28 

U.Zi»_Q 

0.83 ±0.07 



Table 5: Theoretical leptonic widths relatively to the experimental value 
Fe+e- {J 1^) in keV. The experimental numbers are from Particle Data Book. 
The second colum reports the charmonium squared radial radial wave function 
and its derivatives at the origin in GeV3+'. In the widths assign to a given 
known particle, in the computation of the width it was used the experimental 
meson mass rather than the theoretical prediction. 





Vo 


QCD 
(BT) 


Power 
Law 


Log 


Cornell 


Exp. 


J/V(1.5) 
^ l+l- (keV) 












5.55 ±0.14 ±0.02 


l+l- (keV) 


2.84 


2.55 


2.19 


2.00 


2.50 


2.48 ±0.06 


V'(4040) [3^Si\ 
l+l- (keV) 


2.41 


1.84 


1.34 


1.15 


1.77 


0.86 ±0.07 


^(4160) [3^Di\ 
l+l- (keV) 


2.28 


1.73 


1.27 


1.08 


1.68 


0.83 ± 0.07 



Table 6: Charmonium leptonic decay widths in keV. 
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Level 


Final state 


Predicted width 


Measured width 




gg 


12.4 (MeV) 


14 ± 7 (MeV) 




77 


3.26 (keV) 






gg 


9.77 (MeV) 






77 


2.57 (keV) 





Table 7: n^So widths. In the calculation, for 2^5*0 we have used the experimental 
mass, while for the 3^5o we used the predicted mass value. The measured width 
reported is the total particle width. 



Level 


assign. 


Final state 


Predicted width 


Measured width 




Xco(lP) 


gg 


7.22 (MeV) 


10.4 ±0.7 (MeV) (a) 






77 


4.06 (keV) 






Xcl(lP) 


qqg 


86.2 (keV) 


0.89 ± 0.05 (MeV) (a) 




Xc2(lP) 


99 


1.63 (MeV) 


1.55 ±0.11 (MeV) (b) 






77 


921 (eV) 


559 ± 83 (eV) (b) 




hc{lP) 


999 


71.5 (keV) 


< 1 (MeV) (a) 






991 


5.76 (keV) 




23Po 


y(3940) 


gg 


8.27 (MeV) 


87 ±22 ±26 (eV) (a) 






77 


4.66 (keV) 






X(3872) 


qq9 


0.103 (MeV) 


< 2.3 (MeV) (a) 






99 


1.99 (MeV) 








77 


1.12 (keV) 




2iPi 


X(3940) 


999 


84.7 (keV) 








991 


6.82 (keV) 




l^£)i 


*(3770) 


ggg 


0.63 (keV) 


23.0 ± 2.7 (MeV) (a) 






ggg 


0.08 (keV) 






*(3836) 


ggg 


0.27 (keV) 




l^D2 




gg 


0.62 (keV) 





Table 8: P-wave and I?-wave widths. In the calculation for the meson mass we 
used the experimental values of the particle assigned. When no assignment is 
given, for the meson mass we use the theoretical prediction, (a) a total width; 
(b) see ISS] - the value decreases if one includes radiative corrections. 
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